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Abstract 

We study a scenario that a hidden gauge boson constitutes the dominant component of dark matter and 
decays into the standard model particles through a gauge kinetic mixing. Interestingly, gamma rays and 
positrons produced from the decay of hidden gauge boson can explain both the EGRET excess of diffuse 
gamma rays and the HEAT anomaly in the positron fraction. The spectra of the gamma rays and the 
positrons have distinctive features; the absence of line emission of the gamma ray and a sharp peak in the 
positron fraction. Such features may be observed by the FGST and PAMELA satellites. 



1 



I. INTRODUCTION 



The concept of symmetry has been the guiding principle in modern physics. The structure 
of the standard model (SM) is dictated by SU{3)c x SU(2)l x f/(l)y gauge symmetries. The 
electroweak symmetry, SU{2)i x f/(l)y, is spontaneously broken by a non-vanishing vacuum 
expectation value (vev) of the Higgs boson, and massive W and Z bosons are generated. It is quite 
natural in the string landscape that there are many other gauge symmetries as well as discrete ones 
realized in nature, and some of the gauge symmetries may be spontaneously broken, leading to 
massive gauge bosons, as in the SM. 

Suppose that some of the hidden gauge bosons are 'odd' under a Z2 parity, whereas all the SM 
particles are 'even'. If the parity is exact, and if the hidden gauge symmetries are spontaneously 
broken, the lightest parity-odd gauge boson is stable and can be a good candidate for cold dark 
matter. Some recent examples are the T-odd U{\) gauge boson in the little Higgs model with 



T-parity 111 Q, Q, I4] and the KK photon in the minimal universal extra dimension model |15|]. 
However, it is not known yet which of unbroken or (explicitly or spontaneously) broken discrete 
symmetries are more common in the string landscape. If a discrete symmetry breaking is a general 
phenomenon, we may expect that the dark matter is not absolutely stable and ultimately decays 
into the SM particles. 

In this paper we consider a simplest case that there exist a hidden U{\)' gauge symmetry which 
is spontaneously broken and a parity under which only the hidden gauge boson changes its sign. 
We assume that theparity is broken by a kinetic mixing between the U{1)' and the SM U{1)y 
gauge symmetries BaLZD. As a result, the hidden gauge boson decays into the SM particles through 
the parity violating interactions induced by the kinetic mixing. If such a violation is so tiny that the 
lifetime of the hidden gauge boson is much longer than the age of our universe, the hidden gauge 
boson can be the dominant component of dark matter. Furthermore, the subsequent decays of the 
SM particles will form continuous spectra of the gamma rays and the positrons in the high-energy 
cosmic ray. With an appropriate amount of the parity violation, we see that those gamma rays and 
positrons could be the source of the excesses of the gamma ray observed by Energetic Gamma Ray 
Experiment Telescope (EGRET) [8, ^ and the positron flux observed by High Energy Antimatter 
Telescope (HEAT) ilOD, MASS iul and AMS UJ] experiments. 

Recently, the gravitino dark matter scenario was extensively studied in the framework of su- 
persymmetry with i?-parity violation [|l3L Il4fl . The decay of the gravitino into the gamma rays 
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was discussed in Refs. nl3[ 



14ll . Furthermore, if the gravitino is heavier than W boson, it decays 



predominantly into aW or Z boson and a lepton. With the gravitino mass of C(IOO) GeV and 
the lifetime of C(10^^) sec, the gravitino decay can simultaneously explain the EGRET anomaly 
in the extragalactic diffuse gamma ray background and the HEAT excess in the positron fraction 
llSn . Our decaying hidden gauge boson has therefore many parallels with the gravitino 
with _R-parity violation. However, in our model, the decay branching ratio of the hidden gauge 
boson to W boson is highly suppressed, which is significantly different from the gravitino case. 

This paper is organized as follows. In SecUIl we present the effective Lagrangian and Feynman 
rules used in our calculations. In Sec. Hill we calculate the spectra of gamma ray and positron flux 
from the decay of the hidden gauge boson and compare them with the observed data. Sec. HVlis 
devoted to discussions and conclusions. 



II. FRAMEWORK 

We consider a hidden Abelian gauge symmetry f/(l)' and the associated gauge boson B'^, and 
introduce a hidden parity under which B'^ transforms as 

5; - -5;. (1) 

The U{iy symmetry is assumed to be spontaneously broken, so that the gauge boson B'^ has a 
non- vanishing mass, m. 

We assume that the low energy effective theory can be written in terms of the SM particles and 
the hidden gauge boson B'^, and that all the SM particles are neutral under both the 11(1)' and 
the hidden parity. We would like to introduce a tiny parity violating interaction. Among many 
possibilities, we will focus on a kinetic mixing term between U{1)' and f/(l)y, since the kinetic 
mixing has the lowest dimension *^ . The low energy effective Lagrangian can be written as 

C D -\b,,B^^ - \b'^,B'^^ + '-B.^B'^'' + ^m25;5'^ (2) 

where e(^ 1) is the coefficient of the hidden parity violation term, and B^ is the gauge boson 
of t/(l)y, B^^ and B'^^ are the field strength of f/(l)y and respectively: 5^ = - 

*' It is also possible to write down non-renormalizable parity violating interactions suppressed by a large mass scale. 
For the magnitudes of kinetic mixing in our work, however, we can safely neglect such interactions if they are 
suppressed by the Planck scale. 
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d^B^p. In the following discussion we neglect O(e^) terms, since e must be extremely small as 
C(10~^^) for the hidden gauge boson to become dark matter. We will also give a possible origin 
of such tiny kinetic mixing later. 

We can remove the kinetic mixing and bring the kinetic terms into canonical form by redefining 
the gauge fields as [ loj 



B = B-eB 
B' = B'. 



(3) 



Note that we omit the Lorentz index fi hereafter. Taking account of the electroweak symmetry 
breaking, the mass terms of neutral gauge bosons are 



-99' 



-99' 
9'\ 



-gg'e ^ ( \ 



9'^e 



B 



(4) 



\-99'^. 9'^^, 4mVi;2y 

where is the neutral gauge field of SU (2) l, v is the vev of the doublet Higgs field H, g and g' 
are weak and hypercharge coupling constants, respectively. After diagonalizing the mass matrix, 
we can express the gauge eigenstates W^, B and B' in terms of the mass eigenstates Z, A, and A' 
as 



with 



cwZ + swA — cwsw- 



m , 



m 
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B = —swZ + cwA 
B' = sw 



m 
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cw = cos Ow = g/Vg'^ + g 



ml = ty + g'^)v\ 



m 



Sw = sin 9w 

,2 



(6) 



where ttia' and mz are the masses of A' and Z, respectively. One can easily see that A and Z 
are reduced to the ordinary photon and Z boson in the limit of e 0. Since interested values 
of e are extremely small, the kinetic mixing hardly affects any SM predictions of the electroweak 
measurements. 
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The A' interacts with the SM particles through the mixings shown in Eq. (|5]). The Feynman 
rules can be derived in a straightforward way by expanding the SM Lagrangian with respect to e 



02 



2 



2111 . and the results are 



lA'j- 9'^ 2 2 1 2 1 2 



'-kyilb A, lli'7 



z) 



WW A' : - 7"'"^ ^^^^, (7) 

where u, d, e and z/ represent all the up-type quarks, down-type quarks, charged leptons and 
neutrinos, respectively, and g^wz the coupling of WW Z in the SM. 

The hidden gauge boson A' decays into the SM particles through the above interactions. Fig.[T] 
shows the decay branching ratios and the lifetime of A' as a function of rriA'- The down-type 
quark decay modes of A' , A' ^ dd(d = d, s, b), dominate over the other modes for tua' — 100 
GeV, whereas they decrease quickly as tua' increases. For m^/ > 120 GeV, the up-type quark 
decay modes. A' ^ uu{u = u, c), become the largest ones followed by charged lepton decay 
modes. A' t^t~ [i = e, /i,r). This behavior can be understood easily as follows. The A' 
becomes more like the Z boson as m^/ approaches m^, and therefore, the partial decay width of 
A' ff is proportional to gy + g\, where gv(A) is the vector (axial) coupling strength of the 
neutral weak interaction. On the other hand, as m^' becomes much heavier than mz. A' tends to 
not feel the electroweak symmetry breaking, so the branching ratios are insensitive to niA'. For 
"rriA' ^ 350 GeV, the decay mode of A' ti is allowed. The following discussions on the gamma 
rays and the positrons are not significantly modified even for m^' > 350 GeV. 



III. GAMMA RAY AND POSITRON SPECTRA 

The hidden gauge boson A' decays into a pair of SM fermions and W bosons if kinematically 
allowed, as we have seen in the previous section. Due to the subsequent QCD hadronization 
processes, a bunch of hadrons are produced, and in particular, continuum spectra for the photon 
and the positron are formed. Throughout this paper we assume that the A' constitutes the dominant 
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Figure 1: Decay branching ratios and lifetime of A', as a function of tua'- 

component of dark matter. Then the produced photons and positrons may be observed in the high- 
energy cosmic rays. In this section we estimate the fluxes of the gamma rays and the positrons 
produced from the A' decay, and see how they may account for the observed excesses. 



A. Gamma-ray flux 

The gamma-ray energy spectrum is characterized by dN^/dE, the number of photons having 
energy between E and E + dE, produced from the decay of one A gauge boson. The main 
contribution to the continuous spectrum of 7 arises from the 7r° generated in the QCD hadroniza- 
tion process. To estimate the spectrum, we use the PYTHIA [22] Monte Carlo program with the 
branching ratios shown in Fig. [T] We also include the real gamma-ray emission, known as inter- 
nal bresstrahlung, from the charged particles from the A' direct decay whose contributions to 
dN^I dE become important compared to that from 7r° when E ttia' /2. However, the features of 
the gamma-ray flux which we discuss below will not change even though the internal bresstrahlung 
effects are not included. The energy spectra dN^/dE for m^/ = 100 and 300 GeV are shown in 
Fig. [21 It is worth noting that there is no line emission of the gamma rays from the decay of A', 
which is present in the case of the gravitino dark matter. 

There are galactic and extragalactic contributions from the decay of A' to the observed gamma 

I — i 

*^ We are grateful to John Beacom |123|| for bringing up this fact to us. 
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Figure 2: Energy spectra of 7 and e+ generated from the decay of A'. 



ray flux. The flux of the gamma ray from the extragalactic origin is estimated as [|l6L llSQ 



eg 



dy- 



y 



-3/2 



AnmA'TA'Hon]!;' Ji d{yE) A + ^ 



(8) 



where c is the speed of light; and are the density parameters of A', matter (including 

both baryons and dark matter) and the cosmological constant, respectively; pc is the critical den- 
sity; ta' is the lifetime of A'; Hq is the Hubble parameter at the present time; y = 1 + z, where 
z is the redshift, and ycq denotes a value of y at the matter-radiation equality. For the numerical 
results, we use n24\ 



0.1099, fiA//i^ = 0.1326, fiA = 0.742, = 1.0537 x 10~^GeV/cml (9) 



On the other hand, the gamma ray flux from the decay of A' in the Milky Way halo is 



2dJj_ 
dE 



E^ 



dN., 



PhaMdi 



(10) 



^halo ^T^rriA'TA' dE 

where phaio is the density profile of dark matter in the Milky Way, f^^^ Phaio{^)dtj is the average 
of the i nteg ration along the line of sight (los). We adopt the Navarro-Frenk- White (NFW) halo 
profile \m 

= , , w.^l , ,2 (11) 
in our calculation, where r is the distance from the center of Milky Way, = 20 kpc, and po is 
set in such a way that the dark matter density in the solar system satisfies p{rQ) = 0.30 GeV/cm^ 
ll26ll with = 8.5 kpc being the distance from the Sun to the Galactic Center. 
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Figure 3: 7 ray flux predicted from decay of A', shown together with the EGRET data 



In order to compare the EGRET results with the above flux from A' decay, we integrate over 
the whole sky except for the zone of the Galactic plane (i.e. the region with the galactic latitudes 
|6| < 10°). For the background, we use a power-law form adopted in Ref. jlSl] 

^2^1 ^ 5.18 X 10-^E-°-^^9 GeVcm-2sr-isec-\ (12) 

where E is in units of GeV. Fig. [3] shows our numerical results. Recall that there are only two 
parameters in our model, i.e. the mass and lifetime of A', and the peak position only depends on 
the mass while the deviation from background is sensitive to both. The signal of tua' = 100 GeV 
(the red line) peaks at around E = 5 GeV region, which is consistent very well with the observed 
data. For uia' = 300 GeV (the blue line), our prediction can still account for the observed excess, 
although the fit is not so good in some energy region. However, given the large errors in the 
observed data, it may be premature to extract any sensible constraint on the model parameters. 
With more precise data from the FGST experiment, we should be able to have more information 
about the A' mass and lifetime, if the excess is indeed from A' dark matter decay. As mentioned 
before, there is no line emission of the gamma rays in our model, because the production of a pair 



of on-shell 7 from A' is forbidden [|27 



2811 ■ Therefore, there is no secondary peak around the high 



end of signal region, which is a characteristic difference between our prediction and the gravitino 
case. 
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B. Positron fraction 



After being produced from the A' decay, the positron will propagate in the magnetic field 
of the Milky Way. The typical gyroradius is much smaller than the size of the galaxy, and the 
positron will propagate along the magnetic field. However, since the magnetic fields are tangled, 
the motion of the positron can be described by a diffusion equation. Neglecting the convection 
and annihilation in the disk, the steady state solution must satisfy [17] 

V ■ [K{E, f) V /e+] + ^ m, r)U] + Q{E, r) = 0, (13) 

where /e+ is the number density of e"*" per unit kinetic energy, K{E, r) is the diffusion coefficient, 
h{E, r) is the rate of energy loss and Q{E^ r) is the source of producing e"*" from A' decay. In our 
case, 

Q{E^r) P^"^ ^^^^ 



niAiTAi dE 



where dN^+fdE is the energy spectrum of e+ from A' decay obtained by using PYTHIA [1220 
(see Fig. O. We notice that because of the decay channel of A' — > e+e~, there is a sharp peak at 

E = mA'/2. 

The solution of Eq. (fT3l) in the solar system can be expressed as 

feAE) = ^— / dE'G{E,E')^, (14) 
rriA'TA' Jo dE' 



where Emux = rriA' /'2, and G{E, E') is approximately given by lll7h 



1016 



G{E, E') ~ t^e^+^^'' )9{E' - E) sec/cm^ (15) 

where E is in units of GeV, 5 is related to the properties of the interstellar medium and can be 



determined mainly from the ratio of Born to Carbon (B/C) ll29h . We adopt parameters, 5 = 0.55, 



a = —0.9716 and b = —10.012 [|l7h . that are consistent with the B/C value and produce the 
minimum flux of positrons. Finally, the flux of is given by 

(i?) = f /e+ = ^ / dE'G{E, E')-f. (16) 

47r A-RniA'TA' Jo dE' 

In addition to e'^ flux from dark matter decay, there exists a secondary e"*" flux from interactions 
between cosmic rays and nuclei in the interstellar medium. The positron flux is considered to be 
suffered from the solar modulation, especially for the energy below 10 GeV. If the solar modula- 
tion effect is independent of the charge-sign, one can cancel the effect by measuring the positron 



fraction, 

^ '\ ■ (17) 

Indeed, most of experiments measured the fraction of positron flux. To estimate the positron 
fraction, it is necessary to include the flux. We use the approximations of the and e"*" 



background fluxes 130 , 



^ ^ ^ 1 + 11^0-9 + 3.2^2-15 

' ^ ' 1 + 110^1-5 + 600^2.9 ^_ 530^4.2 

^'f(E) = — GeV-^cm-2sec-^sr-\ (18) 

where E is in units of GeV. Therefore, the fraction of e"*" flux is 



(19) 



where is a free parameter which is used to fit the data when no primary source of e"*" flux exists 



320. Note also that the primary flux of e , in the denominator of Eq. (fT9l) should 



include the contributions from dark matter A' decay as well. Our numerical results for niA' = 100 
GeV (the magenta line) and 300 GeV (the green line) are shown in Fig. 31 The prediction of our 
model is consistent with the observed excess quite well, and the position fraction starts increasing 
around ~ 20 GeV and ^ ~ 10 GeV for m^' = 100 GeV and tua' = 300 GeV, respectively. 
Another key feature of the signal prediction is that the curve drops off sharply at E = m^//2 
mainly due to the contribution of e"*" from A' — * e~^e~ decay channel, i.e. the peak seen in Fig.|2l 
These characteristics can be checked by the upcoming PAMELA data. 



IV. DISCUSSION AND CONCLUSIONS 

Let us here discuss briefly how such a small kinetic mixing in Eq. Q can arise. If there is 
an unbroken Z2 parity symmetry under which the hidden gauge boson flips its sign, the kinetic 
mixing would be forbidden. This parity may arise from a more fundamental symmetry or it may 
be an accidental one. In order to have a non- vanishing kinetic mixing, we need to break the 
parity by a small amount and transmit the breaking to the kinetic mixing To this end we 

*^ With the presence of the parity symmetry, a small breaking is natural in the sense of 't Hooft. 
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Figure 4: Fraction of e"*" flux from A' decay, shown togetlier witli experimental data U. 
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introduce messenger fields and ^2 that have both f/(l)y and f/(l)' charges as (1, 1) and (1,-1; 
respectively. Under the parity, they transform as: 



6 ^ 6 



(20) 



Note that the charge assignment is consistent with the parity transformation (see the discussion 
below Eq. (l23l)). We assume that the messenger mass scale Mm is very large, e.g. the grand 
unified theory (GUT) scale. Suppose that the parity is broken in the messenger sector in such a 
way that ^1 and ^2 obtain different masses. Let us denote the mass difference by AA/„, which 
parametrizes the amount of the parity violation. Integrating out the messengers, we are then left 
with a small kinetic mixing, e, given by 

' ~ i^iwr' ^^^^ 

where gh is the coupling constant of the hidden U{1) gauge symmetry. For AM„i ~ C(TeV) and 
Mm ~ 10^^ GeV, we obtain the tiny kinetic mixing of the right magnitude of C(10~^^). 

It is tempting to identify the origin of the parity violation with the spontaneous breaking of 
the U{iy. As an illustration, we will present a toy model below. As we will see later, there are 
some dangerous couplings which may spoil the stability of the A' in this model. However, these 
problems could be solved by embedding the model into a theory with supersymmetry or an extra 
dimension(s). 

Suppose that all the matter fields in the hidden sector are neutral under the parity. Then, the 
gauge field B'^ cannot have any interactions with those hidden matter fields, since they are forbid- 
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den by the parity. Let us introduce two scalar fields, 0i and 02, which transform under the parity 

as 

01 ^ 02, (22) 

and we assume that 0i and 02 are neutral under the SM gauge symmetries. Then the following 
interactions are allowed: 

C D ighB'^{<j)ld^<Pi - h.c.) - ighB'^{4>ld''4>2 - h.c). (23) 

The parity is spontaneously broken if one of the two scalars develops a non- vanishing vev. To this 
end, we consider the following potential: 

A(|0i|2 - vlf + A(|02|2 - vlf + 2«:|0in02r, (24) 

where , A and k are real and positive. For k > \ there are four distinct vacua, (0i, 02) = 
(±f/i, 0) and (0, ±f/i). We take one possibility of them, (0i, 02) = {v^, 0), as an example. There- 
fore, the hidden U{\) gauge symmetry is spontaneously broken by the vev of 0i, and the associated 
gauge boson acquires a mass, m = \/2ghVh, by eating the imaginary component of 0i. 

In order to transmit the parity violation, we introduce couplings between 0i 2 and the messenger 
fields, 

- £ D Mim' + i^r) + ^'m\i? + m'W). (25) 

where Mm is the messenger mass, and k' is a real and positive constant. After the U{1)' is spon- 
taneously broken, masses of and ,^2 are slightly different: m^^ = + k'v^^ and = M^. 
After integrating out these heavy messengers, we obtain the kinetic mixing in Eq. ^ with e given 
by 

_ 9hg' i^'vl 

^ ~ 167r2 Ml ■ 

For gj,^ k' ^ C(l), Vh ~ 1 TeV, and M^ ~ 10^^ GeV, we obtain e ~ 10~26. We can therefore 
realize more or less the correct magnitude of e needed to account for the excesses of the gamma 
rays and the positrons in this toy model. 

In the above toy model, we have assumed that the SM particles couple to the hidden sector only 
through the messenger fields. It is also possible to introduce direct interactions between 0i 2 and 
the visible sector. For instance, we can couple them to the SM fermions as 

C D ^(|0iP + \h\')fLfRH + h.c. , (26) 
12 



where a ~ 0(1), Mp ~ 2.4 x 10^^ GeV is the reduced Planck scale, and Jjk^l) is the right-(left- 
)handed fermion and H is the Higgs field. Although this interaction does not lead to the gauge 
kinetic mixing, it induces the decay of the hidden gauge boson into a fermion pair after breaking 
the parity. However, the decay branching ratio through such interaction is negligible, compared to 
that through the kinetic mixing with e ~ 10^^^. 

More dangerous direct couplings are those between 0i^2 and the Higgs field, ( 1 0i p + 1 02 P) | -f?^ 
If there exist such interactions, the A' could decay into the Higgs bosons immediately. If we extend 
the model into a supersymmetric one, the direct couplings to the Higgs bosons can be suppressed 
with the aid of the _R-symmetry. In a theory with an extra dimension, it is also possible to suppress 
the direct couplings with appropriate configurations of the branes (e.g. the visible particles on one 
brane, while the hidden particles on the other). 

It is also known that many U{1) symmetries appear in the string theory, and the kinetic mixing 
can similarly arise in the low energy theory by integrating out heavy string states that have charges 
of two U{1) gauge symmetries. It has been extensively studied how large the kinetic mixing can be 
in e.g. Ref. i34\ (see also Ref. pSi] and references therein). For instance, in a warped background 
geometry, we can have an exponentially small kinetic mixing [Ssl. 

Let us also have some comments on the the GUT. So far we have assumed the existence of 
the kinetic mixing between U{1)y and the hidden U{iy. If one of the U{1) gauge symmetries 
actually sits within an unbroken non-Abelian gauge symmetry, such kinetic mixing is not allowed. 
This does not necessarily mean that kinetic mixings are incompatible with GUT, because the GUT 
gauge group is spontaneously broken. Indeed, kinetic mixings between the 11(1)' and the GUT 
gauge fields can arise below the GUT scale by picking up non-vanishing vevs of the Higgs bosons 
responsible for the GUT breaking. 

Let us also discuss how the hidden gauge boson could be generated in the early universe to 
account for the observed abundance of dark matter. For simplicity, we neglect a numerical coef- 
ficient of order unity in the following discussion. In the presence of the messenger fields, there 
appears at one-loop level a following interaction between the f/(l)' and U{1)y gauge fields: 

The hidden gauge boson A' will be produced through the above interaction most efficiently at the 
reheating. The A' abundance is roughly estimated to be 



3 ■ 1013 GeV j llO^^Giv' ' ^^^^ 
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where s and Tf{ are the entropy density and the reheating temperature, respectively. For the mass 
niA' ~ 0(100) GeV, a right abundance of A' can be generated from the above interaction for 



Tfj ~ 10 GeV. Also the A' can be non-thermally produced by the inflaton decay [S 
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In the above model, the parity is spontaneously broken by the vev of (pi. If the break- 
ing occurs after inflation, domain walls connecting two of the four vacua in Eq. (I24l) will be 
formed [3^ 



4111 . which can be the cosmological disaster. There are several means to get around 



this problem, and one of which is to introduce a tiny explicit breaking of the parit y sy mmetry. 



Therefore the domain walls are not stable, and eventually annihilate after collisions [|41L 
Another solution is to assume that the breaking occurs before inflation. The last one is to assume 
that the initial positions of the scalars 0i and 02 are deviated from the origin. In the last case, the 
domain walls, if formed, will be annihilated eventually 

In this paper we have considered a possibility that a hidden gauge boson A', which constitutes 
the dominant component of dark matter, decays into the SM particles through the kinetic mixing 
term that breaks the Z2 parity symmetry. As a result, the branching ratios are solely determined by 
the mass of the hidden gauge boson. Continuum spectra of photons and positrons are generated 
from A' — ^ i'^£~ (i = e, /i, r), and from the decays of hadrons, mainly vr"'^, produced in the 
subsequent QCD hadronization process If the mass of A' is about C(IOO) GeV and its lifetime 
is of order (9(10^^) seconds, those gamma rays and positrons from A' decay may account for the 
observed excesses in the extra galactic diffuse gamma ray flux and the positron fraction. Interest- 
ingly, in our model, the spectra of the gamma rays and the positrons have distinctive features: the 
absence of line emission of the gamma ray and a sharp peak in the positron fraction. Such features 
may be observed by the FGST and PAMELA satellites. 
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The significant amount of the gravitational waves may be produced in the collisions of domain walls Il44ll . The 
hidden gauge bosons are also produced by the annihilation processes of the domain walls. 

There also exist antiprotons produced from the A' decay, whose flux can also be calculated and compared with 
the present data. However, we did not pursue the detailed comparison in this paper due to the large experimental 
uncertainties and our poor understanding of diffusion models. 
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Note added: Very recently the PAMELA group reported a steep rise in the positron frac- 
tion 113 311 . which is nicely explained by our model. 
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